Vesicle shape, molecular tilt, and the suppression of necks 
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Can the presence of molecular-tilt order significantly affect the shapes of lipid bilayer membranes, 
particularly membrane shapes with narrow necks? Motivated by the propensity for tilt order and 
the common occurrence of narrow necks in the intermediate stages of biological processes such as 
endocytosis and vesicle trafficking, we examine how tilt order inhibits the formation of necks in the 
equilibrium shapes of vesicles. For vesicles with a spherical topology, point defects in the molecular 
order with a total strength of +2 are required. We study axisymmetric shapes and suppose that 
there is a unit-strength defect at each pole of the vesicle. The model is further simplified by the 
assumption of tilt isotropy: invariance of the energy with respect to rotations of the molecules about 
the local membrane normal. This isotropy condition leads to a minimal coupling of tilt order and 
curvature, giving a high energetic cost to regions with Gaussian curvature and tilt order. Minimizing 
the elastic free energy with constraints of fixed area and fixed enclosed volume determines the allowed 
shapes. Using numerical calculations, we find several branches of solutions and identify them with 
the branches previously known for fluid membranes. We find that tilt order changes the relative 
energy of the branches, suppressing thin necks by making them costly, leading to elongated prolate 
vesicles as a generic family of tilt-ordered membrane shapes. 

PACS numbers: 87.16.Dg, 61.30.Gd, 02.40.Hw, 46.70.Hg 



I. INTRODUCTION 



The interplay of surface curvature and liquid- 
crystalline order finds its fullest expression in the man- 
ifold and complex biological realizations of the bilayer 
membranes surrounding cells and intracellular organelles. 
Helfrich [l| was one of the first to connect membrane 
shape and molecular order, by realizing that the spon- 
taneous curvature of a bilayer membrane could arise 
from the spontaneous splay of the ordered rod-like lipid 
molecules comprising the membrane. But the biolog- 
ical world offers richer varieties of orientational order 
and shape that remain to be understood. For example, 
lipid molecules typically tilt relative to the normal of the 
membrane 0, [3j 13 ; ^^'^ recently become possible 

to image tilt-ordered domains on the surface of curved, 
micron-scale membranes [5]. Furthermore, curvature has 
proven to play an active role in cellular and subcellular 
processes Narrow necks with small mean curvature 
but large negative Gaussian curvature are relevant to bi- 
ological membranes that compartmentalize through bud- 
ding, since this neck geometry allows separate membrane- 
bound compartments to be budded off, while avoiding 
high-energy membrane shapes. Neck formation is uni- 
versal and crucial to the phenomena of endo- and ex- 
ocytosis 0, @], viral entry and budding, the traffic of 
continual fusion and fission of vesicle and Golgi mem- 
brane, and the interconnections between Golgi stacks [9] 
and between the smooth and rough endoplasmic reticu- 
lum. Because of the close association of these phenomena 
with cell function, it is crucial to understand the forces on 
membrane necks and the constraints on their formation. 




FIG. 1: Effect of tilt order on membrane shape. Left: ffuid 
membrane with reduced volume v — 0.706 and spontaneous 
curvature Co = 2.4/7?o, where Ro is an overall length scale 
defined in Sec. IIIII Right: Tilt-ordered membrane with same 
parameters and tilt modulus Km = 2.0k, where k is the bend- 
ing stiffness of the membrane. The arrows represent the tilt 
order. A -1-1 defect sits at both the north and the south poles 
of the vesicle. 



In this paper we numerically calculate the equilibrium 
shapes of axisymmetric vesicles with tilt order. Figure [1] 
shows how dramatic the effect of tilt can be. The two 
vesicles have the same resistance to bending, the same 
enclosed volume, and the same area, but the one on the 
right has tilt order whereas the one on the left has no 
tilt order. As we review below, the tilt order may be de- 
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scribed by a vector field which is tangent to the vesicle 
surface. Molecular interactions prefer uniform tilt order, 
which may be realized on a surface with zero Gaussian 
curvature such as a plane or the surface of a cylinder. But 
it is impossible to have a uniform vector field on a sur- 
face with nonzero Gaussian curvature, such as a sphere or 
the neck connecting the two spheres. Therefore, uniform 
molecular order and Gaussian curvature are incompati- 
ble [l0|, HH ■ In particular, as long as the molecular inter- 
actions are strong enough, the elongated prolate shape 
of Fig. [T] will be preferred over a shape with a neck. 

We begin with a brief discussion of the relation of our 
work to previous work on membrane shapes and orienta- 
tional order. Section [IIII describes our minimal isotropic 
tilt model, coordinates, and numerical method. Our 
analysis and methods are straightforward, but we de- 
scribe them here to make our paper self-contained. In 
Sec. IIVI we present the main results, which are the en- 
ergy as a function of reduced volume for several different 
branches of solutions, and phase diagrams for shapes. In 
the final section we discuss the implications and limita- 
tions of our analysis. 



II. RELATION TO PREVIOUS WORK 

The equilibrium shapes of closed fluid membrane vesi- 
cles have been studied theoretically and experimentally 
for many years (see [H, [3] , and references therein) . In 
the spontaneous-curvature model pioneered by Helfrich, 
a patch of membrane has a resistance to bending but 
is curved in the absence of external loads A differ- 
ent approach known as the bilayer-coupling model ac- 
counts for the bilayer structure of a membrane by im- 
posing a constraint on the number of molecules in both 
monolayers [3, . Both models predict the same set 
of vesicle shapes. However, the spontaneous-curvature 
model predicts that most shape transitions are discon- 
tinuous, while the bilayer-coupling model predicts con- 
tinuous transitions jTB]. To simplify our discussion, we 
will consider the spontaneous-curvature model with tilt 
order. 

As alluded to above, the basic physics governing the 
interaction of vesicle shape and orientational order is the 
incompatibility of Gaussian curvature and uniform or- 
der. This incompatibility is a local property: a patch 
of surface with Gaussian curvature cannot have uniform 
tilt order. The global topology of surfaces also constrains 
the number and strength of point defects in the orienta- 
tion order fleld, via the Poincare-Brouwer theorem, which 
states that the total defect strength of a vector field on 
a surface is equal to the Euler characteristic (see [17] for 
an elementary proof). In our problem, a point defect is 
an isolated point where the tilt order vanishes, and the 
strength of the defect is the number of rotations of the 
tilt order field around that point. For a vesicle with the 
topology of a sphere, the total defect strength is +2. It 
is natural to suppose that the lowest energy states have 



two -1-1 defects at antipodal points (we shall impose this 
two-defect configuration) . 

MacKintosh and Lubensky modeled a vesicle with 
spherical topology made up of molecules undergoing a 
transition from an untilted smectic-^ phase to a tilted 
smectic- C phase [l3| . They found that an initially spheri- 
cal vesicle elongates into a prolate shape, with most of the 
Gaussian curvature concentrated near the defects that 
form at the two poles. They calculated the change in 
shape for this transition assuming fixed area, but they did 
not constrain the enclosed volume. Other work has ex- 
amined the transitions among spherical, cylindrical, and 
toroidal vesicles with orientational order , ag ain without 
the constraint of fixed enclosed volume [liriTsj. In the 
current work, we impose the more realistic double con- 
straint of fixed volume and fixed area, and solve for the 
shape. Also, our numerical method allows us to study 
shapes with large deflections from the spherical geome- 
try. Therefore we can study the effect of tilt order on 
non-spherical shapes such as the pears and oblates pre- 
dicted by the fluid membrane model. Rather than study- 
ing the transition in the tilt order (as in [l3|)) we focus 
on the shape effect: the effect of the tilt modulus Km (the 
elastic constant governing the resistance to non-uniform 
tilt order) on the overall membrane shape. 

Topological defects can also form for geometrical rea- 
sons, even on surfaces such as tori which do not require 
any defects in orientational order. Our work is comple- 
mentary to recent work on the formation and interac- 
tion of such defects on a fixed but arbitrarily curved 
surface [l^ Instead of prescribing the shape and 

solving for the orientational order field, we prescribe the 
positions of two defects and solve for the vesicle shape 
and tilt field. We disallow additional defect formation, 
and discuss the validity and limitations of this restriction 
in Sec. H 



III. THE MODEL AND ITS ANALYSIS 

We make several simplifying assumptions in our anal- 
ysis. Since the vesicles we consider are much larger than 
their constituent molecules, we use continuum mechan- 
ics in the long-wavelength approximation. Thermal fluc- 
tuations are disregarded. We assume that the vesicle 
shapes are surfaces of revolution, and the tilt configura- 
tion is axisymmetric. In particular, the defects required 
by topology are assumed to sit at the two poles of the 
vesicle. These assumptions reduce the partial differen- 
tial equations governing the shape and tilt configuration 
to ordinary differential equations, which greatly simpli- 
fies our calculations. Also, we suppose that the bilayer 
membrane is thin compared to the characteristic size of 
the vesicle. Since the lipid bilayers are approximately 
two nanometers thick this assumption is highly ac- 
curate for vesicles of micron size and larger. A conse- 
quence of this assumption is that stretching is much more 
costly than bending; therefore, we demand that the to- 
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tal area remain constant. And, although membranes are 
permeable to water, osmotic effects resist changes in vol- 
ume leading us to fix the volume. As we explain 
below, we use a minimal model for the orientational or- 
der, disregarding anisotropic couplings between the tilt 
field and the membrane curvature. We also disregard 
chiral interactions, which have been shown to be impor- 
tant in models for lipid tubule structure [2l| [13] and a 
proposed mechanism for budding [H, . 

We will study how vesicle shape depends on area and 
volume in the presence of tilt order. Just as in the case 
of fluid membranes, we will see that the bending energy 
and tilt stiffness energy terms are scale invariant. This 
invariance allows us to vary area A and volume V by 
changing one parameter, the reduced volume v [T3]. The 
reduced volume is the ratio of the actual volume of a 
vesicle to the volume of a sphere with the same area as 
the vesicle. If Rq is the radius of the sphere with area A, 
then V = y/(47ri?g/3). 

A. Parametrization and geometry 

In this section we describe our parametrization and fix 
the notation; see [1^ and [13] for general discussions of 
differential geometry applied to membranes. We choose 
the z axis to be the axis of symmetry of the vesicle and 
represent points on the surface of the vesicle by the three- 
dimensional vector X((/), s) = (r(s) cos </>, r(s) sin 0, z{s)), 
where (f) and r are plane polar coordinates in the xy plane, 
and s is the arclength measured from the north pole of 
the surface along a line of longitude (see Fig. [5]). Define 
tp{s) to be the angle between the tangent vector 9sX 
along a longitude and the horizontal axis. Then dz/dr = 
Zs/ts = — tan'0(s), where < V'(s) < tt. In terms of ip, 
we have Ts = cos'i/' and Zg = — sin^. 

/f Z 




+1 defect 



FIG. 2: Vesicle coordinate system and assumed location of 
defects. The shape is a surface of revolution about the z axis. 
The vector e2 points into the page as indicated by the {^. 

An orthonormal basis in the tangent plane of the sur- 



face is given by 

§1 = = (cos -0 cos (/), cos sin0, — sin?/') 

§2 = S^X/l^^X] = (-sin0, COS0, 0), (1) 

where the s and subscripts denote partial differentia- 
tion with respect to the coordinates s and 0, respectively. 
We construct the outward normal n to the surface using 
the orthonormal frame: 

n = ei X §2 = (sin ip cos (p, sin ip sin 0, cos (2) 

The metric tensor gij of the surface is given by 

g,, = 9.x . a,X =1^1 (3) 

where the indices i,j = 1,2 label the coordinates s and 
(j), respectively. As usual, we denote the inverse of the 
metric tensor by 5'^ , and we use g^^ to raise indices. The 
second fundamental form Kij is defined by 

K.,^n.d.d,X^-^ .sin^)- 

From the second fundamental form we construct the 
mean curvature H and Gaussian curvature K: 

H = -\g'^K,, = 1(0, -I- sin0>/r) (5) 
K = dct{K\j) = det(/i^,j) = i^ssm^/r, (6) 

where repeated indices have been summed over. Note 
from Eq. ([5]) that we use a convention in which the sphere 
has positive mean curvature. 

We assume that the lipid molecules on the surface of 
the vesicle are tilted at a preferred angle with respect 
to the surface normal n, as in a smectic-C phase [3^ . 
Let the vector m denote the projection of the directors 
of the lipid molecules onto the tangent plane. Since we 
do not study the transition between tilted and untilted 
phases, it is convenient to normalize m to make ]m] = 1 
well away from topological defects. In terms of the local 
orthonormal basis of the tangent plane, we have 

m = _B(cos 6*61 + sin 6*62) = i?m, (7) 

where the amplitude B vanishes at defect centers, and 
approaches unity far from defect cores. 

We must assign an energy penalty to nonuniform con- 
figurations of the tilt field m. For a flat surface, m is uni- 
form if the components of m are constant in the standard 
Cartesian basis. Thus, a suitable energy density would 
be proportional to dim^diunK However, on a curved sur- 
face, m can vary with position not only because m* varies 
with position, but also because can vary. Furthermore, 
it is only the tangential component of derivatives of m 
that enter the tilt stiffness terms; normal components 
add to the resistance to bending and therefore may be 
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absorbed in the bending energy term, discussed below. 
These features are captured by the covariant derivative 

Dim = dim — (fi • dim)h 

— {diin^ — m^ili)ei + {diinn? + Tn^^li)e2, (8) 



where the "spin connection,' 



(9) 



is the rate at which the frame {ei, §2} rotates about the 
normal n as the ith coordinate increases. The covariant 
curl of the spin connection is the Gaussian curvature: 



K 



1 

V9 



(10) 



where g is the determinant of the metric tensor and 
is the antisymmetric symbol with £12 = 1 [l7j . We can 
now see why Gaussian curvature is incompatible with 
a uniform tilt field. Consider a tilt field on a surface 
of nonzero Gaussian curvature. A tilt field on a curved 
surface is uniform if it has a vanishing covariant deriva- 
tive. Writing Eq. ([8]) in terms of B and 9, we find that 
Dim = implies diB = and di9 + ^li = 0. Consider 
a patch on a curved surface that includes no defects, so 
that 6 is smooth. Then to solve diO+fli = for 9 we must 
have tijdiVlj = [l3]- Therefore, the tilt field cannot be 
uniform on a patch with Gaussian curvature. 

For our parametrization, = (O,cos^): the frame 
{§1,82} rotates about the normal as changes, but not 
as s changes. Consistent with our assumption that the 
shape of the surface is axisymmetric, we assume that the 
orientation of the molecules is axisymmetric as well, with 
d^B — and S^fl = 0. Therefore, 



D^m — B cos ijjm^ 
D,m = B^m + B9^m 1 , 



(11) 
(12) 



where ihj^ 



sin0ei+cos6'e2, Bs = dgB^ and 9s ~ ds 



B. Free energy 



either monolayer [20|. Since we consider a closed surface 
with fixed topology, the Gauss-Bonnet theorem ensures 
that the integral of the Gaussian curvature is indepen- 
dent of shape and contributes only an overall constant to 
the free energy [33|. Therefore the term proportional to 
KG may safely be disregarded. 

The elastic free energy for the tilt order is a sum 
of many terms, including costs for splay and bend of the 
director field, and many terms coupling the director field 
to the vesicle shape [23, HI]- To simplify our task, we 
demand that the energy be isotropic in tilt, i.e. invariant 
under arbitrary rotations of m about the normal n. This 
symmetry rules out all of the anisotropic terms, leaving 
only a minimal coupling of the tilt order to shape: 



Kr„D,mWr 



+ -(l-TOiTO') ^/gdsc 



(14) 



The first term of Eq. (|14p gives a preference for a uni- 
form tilt field. Since we impose isotropy, the free energy 
density at a point is independent of the direction of m 
relative to the principal directions of curvature. (We dis- 
cuss how anisotropy may affect our results in Sec |V| . 
Note that the Frank elastic constant Km has the same 
dimensions as k, implying that the effects of this term 
are comparable to the effects of the bending term of 
Eq. (fT3|) . The second term of Eq. (|14p gives a prefer- 
ence for |m| — 1. We assume that we are deep in the 
ordered phase, so that XR^/Km ^ 1- The length scale 
-\/ K„i/X determines the radius of the defect core, wherein 
|m| falls steeply to zero. 

The shape of the vesicle and the orientation of the 
tilted molecules on its surface are determined by mini- 
mizing F subject to a given surface area A and volume 
V. To impose these constraints, we introduce Lagrange 
multipliers S and P. It is convenient to treat r and ^ 
as independent variables in the variation of the free en- 
ergy. Therefore, we introduce an additional Lagrange 
multiplier function 7(5) to impose the local constraint 
Ts — cos'i/'- Thus, our task is to minimize 



The free energy F of the vesicle is the sum of terms 
associated with the bending of the vesicle and terms 
associated with the tilt vector order parameter field: 
F — Fh + Fjn- We use the Helfrich model for bending 
energy. 



Fn 



'-{2H - c^f + kgK 



V?dsd0, 



(13) 



In Eq. (I13|) . K is the bending modulus, typically 10- 
Ibk^T and kg is the Gaussian rigidity. The spon- 
taneous curvature cq is twice the preferred value of the 
mean curvature for a patch of membrane. Spontaneous 
curvature of a bilayer membrane can arise either from the 
sum of the inherent spontaneous curvatures of the mono- 
layers, or from a difference in the number of molecules in 



F' = F + + + K / 7(s) (r,. - cos i/')dsd0. (15) 



It is convenient to scale the Lagrange multipliers by k: 
fi = E/k and P = P/k. Then F' can be written as 

F' = 27rK / fi'^,^s,r,rs,B,Bs,9s,j)ds, (16) 
Jo 

where the upper integration limit L is the total arclength 
along a longitude from the north to the south pole, and 

n,p,^s,r,rs,B,Bs,9s,j) 

— (sin^/r + ips ~ Cq)^ + T, + —Prsinip 



5 



2k 



B, 



+^[s){rs — coa-ip). 



(17) 



C. Euler-Lagrange equations 



The Euler-Lagrange equations which extremize F' are 
given by 



Us 

Bs 
Ws 



U sin Tp 



Co, 



— cos V' + 7 sin H — Pr cosih 
r 2 

K^B^ 



KT 



■ cosijjsimp, 



Kmr' 
= -Kb{1-B^) 

K 

KmrB ( cos^ V 



{tBHs)^ 



K 

r' 

2k 

COS?/', 

0. 



(18) 

(19) 
(20) 

(21) 



— cos^ tp 



kW 

KmT 



B 



2q2 



(22) 

(23) 
(24) 



We have introduced two auxihary functions U and W 
to obtain first-order differential equations, which are re- 
quired for our numerical routine. 

The Euler-Lagrange equations can be simplified some- 
what by the observation that 9 is independent of s. 
Equation ([24|) implies that rB^Os is independent of s. 
However, because r and B both vanish at the poles, 
rK'^B^Os = everywhere. Since B and r are nonzero 
only right at the poles, 6's = everywhere. This result 
can be seen more directly by examining the form of the 
free energy F^, Eq. ^1)1 . Since fig = fii = 0, Fm is mini- 
mized when 6 is independent of s. Note that the isotropy 
of the energy means that F is independent of the angle 
0. In Fig. [U for example, we chose to have the directors 
aligned with lines of latitude. The choice was arbitrary 
- directors aligned along lines of longitude, or any other 
direction, would result in the same free energy and same 
shape. 

For the remaining six equations P^H^5|) we have six 
corresponding boundary conditions. The angle ^ and 
radius r are fixed at either pole: 0(0) = r{0) = r{L) — 0, 
and ip{L) = tt. The amplitude B at either pole vanishes 
due to the assumed presence of defects, -B(O) = B{L) = 
0. However, L is still unknown and must be solved for 



along with the shape. To determine L, we reparametrize 
the problem, introducing a new independent variable t 
such that i = at the north pole and t = 1 at the south 
pole: s — Lt, where L is a constant The Euler- 

Lagrange equations (fT8l - [24|) are therefore modified by re- 
placing d/ds with (l/L)d/di. Note that the boundary 
conditions on r, ■0, and B at s = L become boundary 
conditions at i = 1. The additional equation allowing us 
to solve for L is dL/dt = 0. To determine the additional 
boundary condition required by this equation, we con- 
sider the symmetry of the free energy under a constant 
shift in s. The free energy density now takes the form 
F' — 2ttk /dt, where 

Ki^,ipt,r,rt,B,Bt,9t,j,St) 

= stfi^P, i;s/L, r, rjsu B, Bs/suOs/su 7), (25) 

and St = L. Since s{t) does not appear explicitly in the 
free energy density, there is a first integral or Hamilto- 
nian function df /dst which is independent of t. Exam- 
ination of the variation of / with respect to s{t) leads 
to df/dst = at t = 1, implying that the Hamiltonian 
function vanishes everywhere. Writing the Hamiltonian 
explicitly in terms of the dependent variables at t = 
yields the desired boundary condition 7(0) = [l^[2^. 

To complete our specification of the shape equations, 
we describe how we implement the constraints of fixed 
area and volume. It is convenient to regard area as a 
function of t. Define A{t) as the area of the portion of the 
vesicle surface north of the line of latitude corresponding 
to t, with a similar definition for V(t). Then 



At 
Vt 
tt 
Pt 



2iTrL 

TTr^ sin ipL 



0, 



(26) 
(27) 
(28) 
(29) 



where Eqs. ([28)) - l|29|) arise because the Lagrange multi- 
pliers are constant, but unknown and shape-dependent. 
The four boundary conditions corresponding to Eqs. (|26l - 
are A{0) = 0, A{1) = 47ri?§, ^(0) = 0, and 
V{1) — 47ri?Q'i;/3. The scale Rq is set to unity in our 
numerical calculations. Finally, the shape is determined 
by integrating Zt — —Lsintp with boundary condition 
z(0) = 0. 

Although the Euler-Lagrange equations are nonlinear 
and must be solved numerically, it is straightforward to 
determine the form of the amplitude B at the centers 
of the defect cores. For example, near s = 0, we have 
cosV' ~ and r « s. Therefore, equations (^0)1 and (^J) 
reduce to Bessel's equation for s ~ 0: 



s^B, 



sB, 



A 



B = 0, 



(30) 



and thus B on s in the core, which has a size set by 
y^K„i/ X, as mentioned earlier. Similar considerations 
apply to the defect located at the south pole. 
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FIG. 3: Dimensionless free energy F'/k vs. reduced volume v for fluid membrane vesicles with co = for three different values 
of Frank elastic constant: (a) Km ~ 0. (b) Km ~ 0.3k. (c) Km ~ 0.5k (d) Km = n. Solid lines correspond to prolate shapes, 
dashed lines correspond to oblate shapes, and the dash-dotted and dotted lines correspond to stomatocyte shapes. The filled 
circles denote continuous bifurcations. The open circles denote limit points at which the vesicle intersects itself. 



In our numerical approach, we treat the Euler- 
Lagrange equations as a two-point boundary value prob- 
lem, using the MATLAB function bvp4c 30] . We choose 
\/Kra — 1000 to make the defect cores small. To avoid 
the divergences that occur in the equations when r — > 0, 
we solve the equations in the interval 5 < t < \ — 5, where 
5 = 0.001. All boundary conditions are now evaluated at 
t — 5 or t— 1 — 6. We use Taylor series to relate the 
boundary values at these new endpoints to the values at 
t = and t = 1: 

7/;(<5) « i:{0) + MQi)5 = i:tm (31) 

B{5) » B{Q) + Bt{Q)5 = Bt{Q)5 (32) 

ij{l-5) « - ^ - (33) 

B{1~5) » B{l)-Bt{l)5^-Bt{l)5. (34) 

The constants '0t(O), -Bt(O), "04(1), and Bt{l) are un- 
known parameters that may be determined by the nu- 
merical routine using the conditions 

U{6) w L^ljt{S)r{5) + s\n{ilj{5)) + CQr{6) (35) 



U{l-5) w Li!t{l - 5)r{l ~ 5) 

+ sin(V'(l - <5)) cor(l - ^) (36) 

W{5) « LBt{5)K^r{5)/n (37) 

W{1~5) « LBt{l-5)Krar{l-5)/K. (38) 

from Eqns. (HH) and We also have 

7(<5) = 7(0)+7t(0)<^ 

~ l{p)+lt{5)5^^t{5)5, (39) 

with an error of order 5^ . The parameter 74 [5) is given 
in terms of the variables at i = (5 by Eqn. (22). 

IV. RESULTS 

For purposes of comparison with the case of a fluid 
membrane with no tilt order, we have repeated the cal- 
culations of Ref. [1^ for spontaneous curvature cq = 
and coi?o = 2.4. At v = 1, the vesicle shape is always 
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FIG. 4: Dimensionless Gaussian curvature KRq versus the 
dimensionless parameter t for oblate shapes with co = and 
various values of Km/ k,- Note that t = at the north pole, 
t = 1 at the south pole, and t = 0.5 at the equator, which lies 
in the horizontal plane midway between the two defects. 



spherical, with bending energy Sttk. First consider the 
case Co = 0. As w is decreased from unity, the lowest 
energy shape becomes prolate, elongating continuously. 
Upon further decrease of v the shape changes discontinu- 
ously to oblate, and finally to the stomatocyte shape. At 

= where the dash-dotted line in Figure[3^ terminates, 
the stomatocyte shape consists of two concentric spheres, 
spaced infinitesimally close, connected by a vanishingly 
thin neck of zero mean curvature, with a total bend- 
ing energy IGttk. Figure [3^ shows some of these shapes, 
and the free energy vs. v over the range for which the 
transitions occur. Note that there are two stomatocyte 
branches; the upper one (dotted line) is metastable. As 
the reduced volume is decreased along the upper stom- 
atocyte branch, the vesicle height along the z axis de- 
creases and the shape becomes more symmetric about the 
horizontal plane midway between the two defects. The 
two branches join at the filled circle, where the stom- 
atocyte shape becomes oblate. Figures [3b and c show 
the effect of tilt order. As Km/n increases, the prolate 
branch of solutions has lower free energy than the oblate 
branch for a greater range of reduced volume, until even- 
tually the oblate branch becomes completely metastable. 
To understand why tilt order favors prolate shapes over 
oblate, note that as v decreases, the prolate shapes ex- 
tend more along the z axis and become more cylindrical, 
leading to a greater region of small Gaussian curvature 
and approximately uniform tilt order. The oblate shapes 
also have regions of small Gaussian curvature, but these 
regions are confined to narrow bands near t « 0.2 and 
t w 0.8 (Fig.SD. 

Figure ID also shows that our assumption that the de- 
fects are constrained to lie at the poles of the vesicle 
is reasonable. Defects of positive sign prefer regions of 
positive Gaussian curvature [3l| . Figure [3] shows that 



for Km/ K < 0.5, the Gaussian curvature at the equa- 
tor is just greater than the Gaussian curvature at either 
pole. Thus, for Km ^ 0.5k, we expect that the configura- 
tion with two defects at antipodal points on the equator 
might have smaller energy than the configuration with 
defects on either pole. Our simplifying assumption of 
axisymmetry prevents us from investigating this possi- 
bility. However, as Km/ n is increased to unity, the re- 
gions of maximum Gaussian curvature lie at the poles, 
and we expect that the minimum energy configuration is 
to have the defects lie at the poles. For the prolate shapes 
and the pear shapes considered below, the poles are al- 
ways regions of maximum Gaussian curvature. Even for 
the stomatocyte, the north pole is the region of maxi- 
mum Gaussian curvature, and the Gaussian curvature is 
roughly constant over much of the inner surface of the 
pocket. 
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FIG. 5: Phase diagram for lowest energy vesicle shapes for 
Co = 0. The lowest energy shapes are shown for each value of 
V and Kmj 1^- The dashed line denotes the line of limit points 
for the stomatocyte shapes. 

A second qualitative effect of tilt order on the co = 
shapes is that the lower stomatocyte branch develops a 
limit point corresponding to self-intersection at the sites 
of the defects. These limit points are denoted by open 
circles in Fig.O) and c. Figure[5]shows the phase diagram 
for the lowest energy vesicle shapes for the case Cq = as 
a function of the Frank constant. The dashed line shows 
the values of Km/n and v for which the stomatocyte 
shape intersects itself. For values of v below the dashed 
line, the prolate branch again becomes the lowest energy 
branch of solutions. 

Turning to the case cqRq — 2.4, we see from Ref. [l^ 
that the fluid membrane shapes are dominated by prolate 
shapes and pear shapes. Fig. [5^. Note that the energy 
scale on the vertical axis is much less than that of Fig. [3^ 
since the mean curvature required by the constraints of 
fixed volume and area is close to cq . The solid dots again 
denote continuous transitions between pear and prolate 
shapes. The effect of tilt order is to increase the energy 
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V 

FIG. 6: Dimensionless free energy F' / n vs. reduced volume v for fluid membrane vesicles with co-Ro = 2.4 and three different 
values of Frank elastic constant: (a) Km/n = 0, (b) K,n/n — 0.08, (c) Km/f^ = 2.0. Solid lines correspond to prolate shapes, 
dashed lines correspond to asymmetric pear shapes, and the dotted line corresponds to symmetric pear shapes. The filled 
circles mark continuous bifurcations. There is a symmetric pear branch in (a), just to the right of the first bifurcation, but it 
is too short to be seen. 



cost of the pear shapes relative to the prolate shapes, 
eventually making the prolate branch of solutions the 
lowest energy branch for all v. It is interesting to note 
that tilt order does not completely rule out pear shapes 
with narrow necks. The phase diagram of lowest energy 
shapes in Fig. [7] shows that pear shapes are allowed in 
a range of v for sufRciently small Km- The neck of the 
pear shape becomes wider as K„i increases. For these 
shapes, there is a slight reduction in the order i? in a 
narrow band around the neck. 



V. DISCUSSION AND CONCLUSION 

Our detailed, systematic calculations show that tilt or- 
der suppresses necks and favors elongated prolate shapes 
for sufficiently large tilt modulus Km- We calculated the 
free energy as a function of reduced volume for several 
branches of solutions, and showed how increasing the tilt 



modulus increases the energy of the non-prolate branches 
relative to the prolate branches, finally leading to a single 
family of prolate shapes. 

Our calculation was based on several important as- 
sumptions. We ruled out many of the possible terms in 
the free energy by taking the free energy to be invariant 
under arbitrary rotations of rh about the normal n. A 
natural way to remove this assumption and study the ef- 
fects of anisotropy without introducing an unmanageable 
number of terms would be to replace Fm with Fd, the 
one-Frank-constant approximation for the orientational 
free energy of the directors d making up the membrane 
surface l2ll: 



Fd 



Km 

2 



(40) 



To compare Fd with Fm, write d = aN + m, where 
N is the local unit surface normal, and a is determined 
by |d| — 1. Rewriting Fd in terms of intrinsically two- 
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FIG. 7: Phase diagram for co-Ro = 2A. The lowest energy 
shapes are shown for each value of v and Km/f^- 



dimensional quantities, we find that 



+ m'K.jKlm^ - 2a (V.m^) K'} ^dsAct). (41) 

The first two terms of the integrand in (|41|1 are isotropic. 
To see that the term K^Kj is already accounted for in our 
minimal isotropic tilt model, note that a matrix satisfies 
its own characteristic equation: 



K5\ ^ 0. 



(42) 



Using g*-' to take the trace of Eq. 



yields K]Kl 



4H^ — 2K. We expect that the term m'^KijK-jjn!' leads 
to a preference for the tilt to order along the long axis 
of a prolate shape, and an increased bending stiffness 
along the direction parallel to rh. The last term in the 
integrand of (|4T|) is like a spontaneous curvature given 
by a gradient of the projected director. Except for the 
fact that the anisotropic terms will lead to a preferred 
direction of rh relative to the principal directions of cur- 
vature, we expect that including these terms would not 
qualitatively change our results. A more dramatic change 
is expected if chiral terms like m^eijKlm'' are allowed. 
This term can give the vesicle a chiral shape, which is 



necessarily non-axisymmetric. It would be interesting 
to add this term alone to our minimal model and study 
vesicles shap es, as has recently been done for tubules and 
ribbons [Ij. 

The most severe assumption of our model is the as- 
sumption of fixed defects at the poles. We have already 
mentioned that positive sign defects prefer regions of pos- 
itive curvature, which suggests there may be lower en- 
ergy, non-axisymmetric shapes than the ones we consider 
here. Furthermore, recent calculations have shown how 
a pair of defects of opposite sign can be pulled apart 
whenever there is a change in sign in the Gaussian cur- 
vature, with the positive defect migrating to the region of 
positive Gaussian curvature, and the negative defect mi- 
grating to the region of negative Gaussian curvature [20[ . 
The magnitude of the curvature in both regions is im- 
portant for determining whether or not a defect pair will 
unbind. Although we completely disregard this effect, we 
expect it will not play much of a role in the low-energy 
prolate shapes, since for those shapes there are at most 
only narrow bands of mildly negative Gaussian curvature 
where either end starts to bow out. An important exten- 
sion of our calculation would be to lift the assumption 
of axisymmetry, and allow both the defect position and 
number and vesicle shape to vary. It would be interesting 
if the techniques of Ref. [s^ could be used to experimen- 
tally determine defect position on vesicles with tilt order. 

Finally, we have systematically studied the effect of 
changing the tilt modulus Km while deep in the ordered 
phase. While we have not presented results on the change 
in shape due to a phase transition from an ordered tilt 
phase to a disordered fluid phase, as in Ref. [Tc|, we ex- 
pect that the effect of increasing the order should be qual- 
itatively similar to increasing Km- 
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